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Mixing is one of the basic functions which automated lab-on-chips require for the effective man- 
agement of liquid samples. In this paper we report on the working principle, design, fabrication and 
experimental characterization of a soft-lithographed micromixer for microfluidic applications. The 
device effectively mixes two liquids by means of chaotic advection obtained as an implementation of 
a Linked Twisted Map (LTM) . In this sense it is chaotic. The liquids are electrokinetically displaced 
by generating rolls through AC electroosmosis on co-planar electrodes. The device performance has 
been tested on dyed Dl-water for several voltages, frequencies and flow-rates, displaying good mixing 
properties in the range of 10 4- 100kHz, at low peak-to-peak voltages (~ 15 4 20 volts). Low voltage 
supply, small dimensions and possibility of fabrication via standard lithographic techniques make 
the device highly integrable in lab-on-a-chip platforms. 
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I. INTRODUCTION 

Since the beginning of the 1990s many research efforts have been devoted to the development of integrated miniature 
fluidic systems for the automated handling of biosamples on chip (lab-on-chip and micro total analysis systems 
(fiTAS)). In many biological and chemical experimental procedures the ability to control the mixing index between 
different reactants is often of paramount importance for fast and efficient reactions. The development of micromixers 
for microfluidic applications has attracted the attention of many research groups worldwide, and a quite large number 
of papers already appeared which report on the development of micromixers based on different actuation strategies 
(see for instance pQ). 

Mixing occurs when at least two species molecules spread around and diffuse uniformly in the bulk of a medium. 
Diffusion is a spontaneous process: the simplest micromixer is a simple channel into which two liquids are injected. 
If D is the diffusion coefficient of a certain species, w the characteristic dimension of the channel cross section, and u 

the average flow speed, the Peclet number is Pe — j)U and the given species mixes homogenously only after having 
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travelled a distance over the channel I ~ jj u - I n order to reduce mixing times, channels with small cross sections 
are required. Nonetheless, the diffusive mixing length may be quite long. For example, mixing large macromolecular 
species, such as DNA or proteins, which typically have a diffusion coefficient on the order of I0~ 11 m 2 s~ 1 in a 200/im 
wide channel, with an average flow speed of 1mm/ s, requires a channel length of about lm with a corresponding 
mixing time exceeding Ibmin. This example shows that passive micromixers, which rely only on spontaneous diffusion 
of species suspended in liquid carriers, are not really compatible with miniature chips. On the contrary, active 
micromixers are able to both shortening mixing length and reducing mixing time by "stirring" the liquids, so that 
their interfacial area increases while the actual diffusion distance is greatly shortened. 

Since flows in microchanncls are typically controlled by viscosity ("creeping flows" at low Re numbers), stirring 
cannot be achieved by means of turbulence. On the contrary, convective transport can be effectively used to generate 
chaotic advection, that is streamlines stretching-and-folding in such a way that originally close fluid particles follow 
trajectories which diverge exponentially in time. 

Also chaotic micromixers can be passive or active. In the former class, chaotic advection may be induced by means 
of geometrical features of the mixing channel ([3] and [5]). In the latter class, active mixing is generated by means 
of an external actuation which can be mechanical (e.g. in [1] self-assembled magnetic microstirrers are used), or 
electrokinetic (e.g. [H |6j [7]). The second type is of greater practical interest because it does not involve moving 
parts. On the other hand, current electrokinetic devices show a few drawbacks due to the typically high electric 
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voltages required. The device presented here overcomes such limitations by taking advantage of low-voltage AC 
electrokinetic rolls generation (see below). 

The device was designed according to the Linked Twisted Map (LTM) strategy described in [5] and based on 
dynamic system theory. As the liquids flow downstream, they encounter a number of mixing elements which induce 
transverse circulating flows. As a fluid particle transverse one element's length, its coordinates in the cross-section 
plane change to new values according to a map of the cross section into itself. The flow through many identical 
elements, then, can be modelled as the repeated application of that map. The mixer would perform well only if the 
map has effective mixing properties in the sense of dynamic system theory. The LTM is one of the possible mixing 
maps which might be implemented in real microfluidic mixers. The main idea of the LTM is that of alternating 
an asymmetric (with respect to one of the channel symmetry planes) circulating pattern with its symmetrical one 
(with respect to the same plane) [8]. The problem has been studied mathematically jS], to give the conditions on 
the circulating flows' "strength" and "overlap", which ensure the map to be mixing. The device proposed in [3J, 
is a striking example of chaotic advection induced by an LTM. In that device, also named staggered herringbone 
micromixer the circulating flows are generated by means of opportunely defined grooves on the top of the channel 
and, as such, it is a passive device. 

Recently [TS] , Sasaki et al. [TU] have demonstrated that electrokinetic rolls can be used for developing an active 
micromixer. In this paper we further develop the approach presented in |10j in the framework of the chaotic LTM- 
mechanism and we present a novel device which differs from that of [10 both in geometry and fabrication technology. 
In particular, we show that the device can be effectively fabricated using glass an polymers, instead of silicon, using 
a low-cost rapid prototyping techniques called soft-lithography [TT] . Electrodes have been fabricated by sandwiching 
Ti-Au-Ti layers, the first one (Ti) serving as an adhesion layer over glass, the second one (Au) to reduce the electrical 
impedance and the third one (Ti) as a dielectric layer, since Ti spontaneously form a coherent dielectric oxide layer. 
Compared to the device described in [10] . current mixer performed well with very high Peclet numbers in the range 
0.6 * 10 6 - 1.8 * 10 6 . 

In Sec. Illlwe report on the design of the micromixer, whose fabrication is described in Sec. Ill The last paragraphs 



are devoted to described the experimental measurements carried on the device (Sec. IV I and to discussion and 
conclusion (Sec. W\. 



II. DESIGN 



The generation of rolls by means of AC electroosmosis has been thoroughly studied both experimentally and 
theoretically |121 1131 1141 115] and it is becoming increasingly popular in the microfluidics community. Examples of 
the employment of electrokinetic rolls include micropumps [16 j and control of microparticles [17] • The mechanism of 
electrokinetic roll generation is schematically illustrated in Fig. [T] Two coplanar micro-electrodes, spaced by a thin 
gap (on the order of 10 fim wide) are addressed with an AC low voltage (< 20 volts). The electric field interacts with 
a ionic diffuse layer, which forms on top of the electrodes, giving rise to a constant slip velocity which points from the 
electrode gap outwards (see Fig. [I] top). Since the liquid is confined in a closed cavity, the slip- velocity results in the 
formation of a couple of counter-rotating rolling patterns (see Fig. [T| bottom). We used asymmetric microelectrodes 
to produce a sequence of asymmetric rolls in order to implement a LTM. The device comprises a T-channel with two 
inlets for the fluids to be mixed and one outlet for the mixed flow. The flow is pressure driven. As the two liquids 
enter the main channel their flow is laminar (Re < 5) with a very slow spontaneous mixing rate. The mixing tract, 
hosting the 2 micro-electrodes on its floor, starts at entrance of the main channel. The electrodes are patterned so as 
to define 5 mixing elements (see Fig [2]), each comprising two zones. In the first zone the gap (20^tm wide) between 
the electrodes is placed at a distance w/3 from one channel side and 2w/3 from the other side, where w — 200^tm is 
the channel width. In the second zone the gap is located symmetrically with respect to the channel axis, so that the 
positions of the two electrodes are interchanged (twisted). Therefore the electrode gap has the shape of a square wave 
moving along the channel axis. Each mixing element generates two couples of asymmetrical rolls, one for each zone. 
The rolls generated in the second zone are identical to those generated in the first zone, except for their interchanged 
(twisted) location. Such spatial twisting is at the basis of the implemented LTM strategy. 

Because of the very low Re number, the resulting flow in the channel can be imagined as the supcrimposition of 
an axial pressure driven flow with transversal and twisted circulating flows, as schematically represented in Fig. [3] 
The device might also work with other periodic shape of electrodes. As an example, in |10j a device with a sinusoidal 
gap between electrode is presented. 
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ELECTRODES ■ 
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FIG. 1: Schematic of electric field lines, slip and tangential velocities (top) and rolls (bottom) on top of planar electrodes 
addressed by an AC voltage. 

observation 




electrode 

FIG. 2: Microphotograph of the micromixer. Black solid lines have been superimposed for a better visualization of the T- 
channel. 

III. FABRICATION 

The channel has been fabricated in PDMS using soft-lithographic techniques [TT]. The electrodes where patterned 
on a glass substrate using a lift-off process. The electrodes were made of an Au layer (~ lOOnm) sandwiched between 
two Ti layers (~ lOnm bottom, ~ 20nm top) as described in [T^]. The bottom layer serves as an adhesion layer, 
whereas the top layer serves for insulation purposes (Titanium develops an insulating superficial thin oxide layer at 
room temperature and pressure conditions). The PDMS channel was aligned and mounted on the patterned glass 
substrate to define the micro-device. Sealing between PDMS and glass was obtained by assembly the parts in a 
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Pressure driven 
flow 



FIG. 3: Schematic representation of circulating transverse flows (rolls) generated by AC electroosmosis which superimpose to 
the axial pressure driven flow within the micromixer. The positions of the asymmetric vortexes are twisted along the channel 
length in such a way as to implement a linked-twisted-map. 

purposely developed mechanical support. The main channel is 200/im wide and 50fim high. The mixing region is 
5mm long for a total of five 1mm long mixing elements. 

IV. EXPERIMENTAL 

The micro-device was treated in an oxygen plasma (3mm, 18W) to make the channels hydrophilic, then it was 
primed with DI water. The two inlet reservoirs were respectively loaded with DI water and methylene blue dyed water. 
The pressure was applied by connecting the inlets to a pressurised air chamber. The pressure was measured using a 
digital manometer (model 840081 by Sper Scientific). The pressure in the air chamber was generated hydrostatically 
connecting the chamber to a vertical column filled with water. The device outlet was left at room pressure (latm). 

The electrodes were connected to a wave form generator with 20V p - p maximum output voltage and addressed with 
a sinusoidal wave. In the experimental tests we observed the region of the channel which follows the mixing region 
using a CCD microscope (from Hirox inc.). Digital videos at 30 frames per second were stored in a PC equipped with 
a frame-grabber (VCE-pro by Imperx Inc.) for off-line analysis. 

Figure [4] shows two pictures (300 x 550 pixels) taken at constant applied relative pressure (0.25psi, corresponding 
to Pe — 0.58 x 10 6 ) when the device is off (top) and on. In both cases the applied voltage is 20V PP at lOOfcf/z. The 
picture shows how the blue dye spreads homogeneously through the channel cross section. 

To quantify the mixing level starting from images analysis we adopted the standard deviation method as described 
in [5] . The frames were first converted to grey-scale and then pixel values / were normalized to 1 , so that corresponds 
to white and to black. Finally, standard deviation a =< I 2 — < I > 2 >!/ 2 of pixel values was computed. 

Two theoretical limit cases exist. For perfectly unmixed fluids [20J half pixels are black and half are white. In this 
case it is a = 0.5. For fully mixed fluids all the pixels have the same grey color and it is a = 0. Based on this 
observation, we define the mixing index as: 

m = 1 - 2a, (1) 

in such a way that m = corresponds to completely unmixed fluids and m = 1 to completely mixed ones. It is worth 
underlining that this definition of mixing index contains an element of arbitrariness as it depends on the contrast of 
the image. Further, since the standard deviation is used, m is affected by the CCD noise. As a consequence, images 
need to be pre-processed before computation of the standard deviation. We proceeded as follow. Since the mixing 
index m is evaluated over the cross sections of the channel, which correspond to lines perpendicular to the axis as 



FIG. 4: Normalized pictures of unmixed streams (top) and mixed streams (bottom) in the micro-channel. The images show 
the region of the channel right downstream with respect to the mixing region. The device was turned off when the top image 
was shot and was turned on at 100kHz, 20V P - P in the bottom image. The Peclet number was 0.58 x 10 6 . 

viewed in the grabbed frames, pixels intensities have been averaged along rows, so as to reduce each 300 x 550 image 
matrix to a 300 x 1 array. This array contains the mean value of pixel intensities along each row. To overcome the 
arbitrariness due to contrast we took the first frame acquired during the experimental session, reduced it to an array as 
explained above, evaluated its maximum and minimum pixel values I m i n , I ma x an d adjusted the contrast by linearly 
mapping the pixel value interval [I m in> Imax] into the interval [0,1]. The same transformation was then applied to all 
the 300 x 1 arrays obtained from the frames acquired during the experimental session. Finally the mixing index was 
calculated for each array. 

We performed a number of experimental tests at different applied pressure differences (AP = 0.25,0.50,0.75 psi, 
corresponding respectively to measured Peclet numbers of Pe — 0.58 * 10 6 , 1.22 * 10 6 , 1.86 * 10 6 ), for different values 
of driving frequency (ranging from 10 2 to 10 6 Hz) and four values of peak-to-peak voltage amplitude (5,10, 15 and 
20 volts). The chosen Peclet numbers correspond to average flow velocities of 29, 61 and 93 mm /sec and flow rates 
of 0.29, 0.61,0. 93/LtZ/s. 

The tests consisted in turning the device on and off while recording a video. Fig. [5] shows a typical temporal 
evolution of the mixing index. For each test we computed the temporal evolution of the mixing index. The average 
mixing index calculated over the high plateau has been taken as a measure of the mixing performance of the device 
at the corresponding values of frequency, voltage and flow-rate. The error on this measure is given by the standard 
deviation of the mixing index on the plateau. 

V. RESULTS AND DISCUSSION 

The results of all tests are collected in Fig. [6] The mixing index is higher for higher applied voltages. All the curves 
display a similar behavior: the mixing index grows with frequency up to the value of 10 kHz, where a plateau begins. 
In most cases a decrease is observed around lMhz, where electrothermal phenomena may occur. Peak performance 
(m ~ 0.85) is observed at frequencies between 10 kHz and 1 MHz, for the lowest Peclet number employed (0.58* 10 6 ). 
The presence of a peak frequency around 100 kHz was expected on the basis of the theoretical analysis of Ref. |18j : 
at lower frequencies the field quickly relaxes into the bulk, whereas at higher frequencies the electric charges do not 
have time to build up over electrodes. 

Measurements also show that the mixing index decreases with increasing flow rates, so that longer mixing lengths 
are required for higher flow velocities. This behavior is expectedly typical for such devices and a similar one was also 
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FIG. 5: Typical temporal evolution of mixing index as the device is turned on and off. 

reported in [TU] . 

We conclude by observing a few advantages and drawbacks of this device with respect to other similar ones reported 
in literature. For what concerns the drawbacks we should notice that, unlike the staggered herringbone micromixer 
[3J, the axial and transverse flows are independent from each-other. Therefore, as the axial velocity increases, one has 
to increase the velocity of the circulating flows accordingly in order to guarantee that the resulting map be mixing. 
Practically this poses a limit to the range of flow rates over which the device performs well. Nonetheless, the devised 
application as micromixer for Lab-on-Chip is not hindered by such limit, since the limit flow-rates are above those 
which can practically be produced in Lab-on-a-chip devices. 

Among the advantages we should mention that unlike other electrokinetic micro-mixers jH |7] which use high 
voltages (100 -i- 1000 volts), here the voltage employed is quite low (< 20 volts). Furthermore, the high frequency AC 
operation preserve the device from observable electrolytic bubbles generation. 

Another nice feature of this device is its scalability. As the channel cross section gets smaller it is expected that the 
same performances be reached with lower voltages. The main advantage of active micromixing is not only that the 
liquid can be stirred at will, but it can be mixed up to a desired level by tuning the applied voltage. This makes the 
device amenable to be integrated with a programmable control system which records pixel values with a CCD sensor, 
process them in real time and tunes the voltage output in order to control the mixing level. 

All these properties, along with the small dimensions and the possibility, which we have demonstrated, of using soft- 
lithography as fabrication technique, makes the device ideal for integration on low-cost, multi-functional Lab-on-a-chip 
platforms. 
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FIG. 6: Mixing performance as a function of driving frequency for different peak-to-peak voltage amplitudes at Pe — 0.58* 10 6 
(A), Pe = 1.22* 10 6 (B), Pe = 1.86* 10 6 (C). Data corresponding to a driving voltage of 100Hz at V pp ^ 15V are not reported, 
being affected by bubbles generated by water electrolysis. 
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